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We study the mixed valence regime of a generalized Anderson impurity model using the bosonization 
approach. This single impurity problem is defined by the U = oo Anderson model with an additional 
density-density interaction, as well as an explicit exchange interaction, between the impurity and 
conduction electrons. We find three points in the interaction parameter space at which all the corre- 
lation functions can be calculated explicitly. These points represent the mixed valence counterparts 
of the usual Toulouse point for the Kondo problem, and are appropriately named the Toulouse points 
of the mixed valence problem. Two of the Toulouse points exhibit the strong coupling, Fermi liquid 
behavior. The third one shows spin-charge separation; here, the spin-spin correlation functions are 
Fermi-liquid-like, the charge-charge correlation functions and the single particle Green function have 
non-Fermi-liquid behaviors, and a pairing correlation function is enhanced compared to the Fermi 
liquid case. This third Toulouse point describes the novel intermediate mixed valence phase we have 
previously identified. In deriving these results, we emphasize the importance of keeping track of the 
anticommutation relation between the fermion fields when the bosonization method is applied to 
quantum impurity problems. 



PACS numbers: 71.27. +a, 71.10+x, 71.28.+d, 74.20.Mn 



I. INTRODUCTION 

The mixed valence problem is a classic problem in con- 
densed matter theory. It describes an impurity with 
three configurations, having both charge and spin de- 
grees of freedom, coupled to a conduction electron bath. 
It differs from the Kondo problem in that low lying local 
charge fluctuations coexist with local spin fluctuations. 
The Kondo problem has been studied by a variety of 
techniques, and is by now well understood. As long as 
the effective Kondo coupling is antiferromagnetic, the low 
energy behavior is described by a strong coupling, Fermi 
liquid fixed point The mixed valence problem 

has also been studied extensively. A variational study by 
Varma and Yafet J5j, and renormalization group stud- 
ies of Haldane j(| and Krishnamurthy et al. , have all 
found that the low energy behavior of the mixed valence 
problem is described by a strong coupling, Fermi liquid 
fixed point. This fixed point is qualitatively similar to 
that of the Kondo problem, though quantities such as the 
Wilson ratio are modified. These works were followed by 
extensive studies on the proper description of the result- 
ing Fermi liquid states of both the Anderson impurity 
and Anderson lattice problems. These later works used 



techniques ranging from Gutzwiller variational wavefunc- 
tions to the slave boson large N expansion [|| . They have 
provided qualitative, and sometimes quantitative, under- 
standings of the physical properties of the heavy fermion 
metals. 

In a series of papers we have revisited the mixed va- 
lence problem in an attempt to identify metallic non- 
Fermi li quid p hases in a two band extended Hubbard 
model P ^O|Jll| . In the limit of infinite dimensions, this 
extended Hubbard model can be solved through a gen- 
eralized Anderson impurity model with a self-consistent 
electron bath. This single impurity model is defined by 
the U = oo Anderson model with an additional density- 
density interaction, as well as an explicit exchange inter- 
action, between the impurity and conduction electrons. 
Away from half-filling, the self-consistency condition im- 
plies that, the associated impurity problem is in the 
mixed valence regime over an extended range of densi- 
ties. The persistence of the mixed valence behavior for 
the self-consistent impurity problem associated with a 
lattice model in infinite dimensions is in fact quite gen- 
eral and is not restricted to the two band extended Hub- 
bard model. For example, in the one band case, the 
average occupation of the impurity in the self-consistent 
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Anderson model is simply the lattice density, which is not 
very close to one (local moment regime) or zero (empty 
orbital regime) over a wide range of densities. This, to- 
gether with our finding that away from half-filling the 
self-consistent conduction electron bath has a finite den- 
sity of states at the Fermi level, imply that the classifi- 
cation of the fixed points in the mixed valence regime of 
the generalized Anderson impurity model with a regular 
electron bath also specifies the possible metallic phases of 
the extended Hubbard model in infinite dimensions. Of 
course, in general, the couplings of the effective Anderson 
model at the cutoff energy scale of the universal regime 
can be different from the bare atomic interactions, just 
like the Coulomb pseudopotential at the Debye tempera- 
ture scale is different from the bare Coulomb interactions 

We studied the generalized Anderson impurity model 
by extending Haldane's renormalization group scheme 
such that the local charge fluctuations and local spin fluc- 
tuations are treated on an equal footing. In the mixed 
valence regime, we found three, and only three, kinds 
of fixed points, which we termed the strong coupling, 
weak coupling, and intermediate phase, respectively. The 
strong coupling and weak coupling phases are the direct 
analogs of their counterparts in the Kondo problem. The 
intermediate phase is entirely new and occurs only in the 
mixed valence regime. Its existence came as a surprise. 
In this new phase, spin and charge excitations are sep- 
arated. From the renormalization group analysis, it is 
expected that the spin-spin correlation functions remain 
to have the Fermi liquid form, while the charge-charge 
correlation functions and single particle Green function 
have an algebraic behavior with interaction-dependent 
exponents. 

The renormalization group procedure is based on a 
Coulomb gas representation for the single impurity prob- 
lem. As usual, the Coulomb gas analysis uses a dilute 
instanton expansion. In the strong coupling and interme- 
diate phases, some of the fugacities flow towards strong 
coupling, and the behaviors of the correlation functions 
can only be inferred through "extrapolating" the scaling 
trajectories beyond the dilute instanton regime. For the 
intermediate phase, this procedure does not allow an ex- 
plicit determination of the exponents. Finally, it is in 
principle possible that additional fixed points, not cap- 
tured by the dilute instanton expansion, may occur. An 
example for the latter arises in the related, though qual- 
itatively different, problem of tunneling through a point 
contact in a Luttinger liquid |l3| , |l4| ]. 

To address these issues, here we study the mixed va- 
lence regime of the generalized Anderson model at some 
particular values of the interactions where the model is 
exactly soluble (in the sense that will be made precise 
below flql ). These points in the interaction parameter 
space are the mixed valence counterparts of the usual 
Toulouse point of the Kondo problem |I6p[| , and are nat- 
urally called the Toulouse points of the mixed valence 
problem. We identify all the possible Toulouse points us- 



ing the bosonization approach. We construct an effective 
Hamiltonian, and determine the single particle, spin-spin, 
and charge-charge correlation functions, at each of these 
Toulouse points. There are two subtle aspects associ- 
ated with these Toulouse points. First, at the Toulouse 
points some of the interactions are larger than the con- 
duction electron bandwidth. This makes the interpreta- 
tion of these Toulouse points subtle. We clarify this issue 
through a comparison with an atomic expansion of the 
original impurity Hamiltonian. Second, when applying 
the bosonization method to the mixed valence problem, 
it turns out to be essential to keep track of the anti- 
commutation relation between the fermion fields. When 
properly understood, the solutions to the generalized An- 
derson model at these Toulouse points substantiate, and 
provide a simple physical picture for, the mixed valent 
phases we have identified through the renormalization 
group analysis. 

Parallel to our previous works, an impurity problem 
defined by the U = oo Anderson model with additional 
screening channels have been studied using the numerical 
renormalization group method |X^ ]. Recently, two other 
groups Jl8|,[l9| have studied an exactly soluble point of 
that impurity model. They have reached conclusions very 
different from ours, and we will comment on the origin 
of these differences. 

The setup of this paper is the following. In Sec. O, 
we review the construction of the Toulouse points in a 
context where the physics is well understood, the Kondo 



model. In Sec. Ill, we summarize the main results of 
our Coulomb gas analysis of the mixed valence regime 
of the generalized Anderson model, and set up the for- 
malism for our Toulouse-point analysis. The Toulouse 
points are discussed in detail in the following sections. 
The Toulouse points discussed in Sec. IV and Se c. M lie 
deep in the strong coupling limit. That of Sec. |VI| de- 
scribes the non-Fermi liquid, interme diate phase. There 
is one more point, presented in Sec. VII, which has an 



effective Hamiltonian corresponding to a rank-2 general- 
ization of the Emery-Kivelson resonant level model pc| ] 
associated with the two-channel Kondo problem. This 
effective Hamiltonian is not exactly soluble. We con- 
clude with a comparison of our results with those of 
related works, and a discussion about realizing the in- 
termediate phase in other models of strongly correlated 
electron systems. Two appendixes are included. Ap- 
pendix A summarizes the bosonization procedure rele- 
vant to our discussion, with an emphasis on the Klein 
factors that keep track of the anticommutation relation 
between fermions of different spins. Appendix B substan- 
tiates the bosonization results presented in the main text 
with those of an atomic analysis of the original Hamilto- 
nian. 
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II. TOULOUSE POINTS OF THE KONDO 
MODEL 

In this section, we introduce the notation, and re- 
view the bosonization method, in the context where the 
physics is well understood, the anisotropic one-channel 
spin— | Kondo model. We show that there are two val- 
ues of J z , the longitudinal component of the Kondo ex- 
change interaction, at which the problem is exactly sol- 
uble |15|] . One of these is the well-known Toulouse point 
JTol , for which the effective Hamiltonian is the resonant- 
level model |2^| . The other occurs at the infinite value for 
J z . The effective Hamiltonian for this second Toulouse 
point is the spin-boson Hamiltonian [^2 23 close to a 
vanishing spin-boson coupling. 

The anisotropic Kondo problem describes a bath of 
spin— i conduction electrons coupled to an impurity 
which can fluctuate between two states, | |> and | |>. 
The Hamiltonian is 



(i) 



Here, c\. a describes a free conduction electron bath with 



energy dispersion E k , and 4(f) = y—jEt^'^L- 
The impurity is located at r = 0, and is locally coupled 
to the conduction electrons through an exchange interac- 
tion. Jj_ and J z are the transverse and longitudinal com- 
ponents of this exchange interaction, respectively, and 
are allowed to take different values. In Eq. (|l]), we have 
written the impurity spin operators in terms of the Hub- 
bard operators, X a p — \a >< /3\, where a, f3 describes 
the two impurity configurations, | f> an d | j>- The 
following constraint 



A TT + Xu — I 



(2) 



supplements Eq. (hi) 



A. Bosonization 

Given that the interaction occurs at r = only, we 
need to keep only the S— wave component of the conduc- 
tion electrons. This S— wave component is defined on the 
radial axis, r E [0, +oo), and can be further decomposed 
into an outgoing and an incoming components. In a stan- 
dard fashion, we extend to the full axis, x £ (— oo, +oo), 
by retaining only one chiral component, which we denote 
by ip a (x). We can then introduce a boson representation 
for the i/) a (x) field |p5[ . The details of this procedure is 
given in Appendix A. At the origin, 



^(0) - Fl- 



(3) 



Here, a is a cutoff scale which can be taken as a lattice 
spacing. $ CT is the shorthand notation for Q a (x = 0). 
An important point is that, $ CT depends only on the 
q 7^ components of the Tomonaga bosons, b qa and b qa . 
The operator FJ, and its adjoint F CT , are the so-called 
Klein factors. They should be thought of as acting on 
the q — sector of the Hilbert space for the Tomonaga 
bosons. More precisely, they can be defined as the rais- 
ing and lowering operators, respectively, in such a Hilbert 



space 



and b qa for q and, hence, 



These operators are unitary, and anti- 
commute among the different spin species. Furthermore 
they commute with b qa 
also with $cr. 

The Kondo Hamiltonian can be rewritten in the 
bosonized form, H = H a + H± + Hu , with, 

Air J ax ax 



Hi 



Hu 



J± 
Ana 
V2S 



(X n FlF^- 1 *^ + H.c. 



(Xf-f - Xu)(— ^) x=0 — 

irp a ax Zir 



(4) 



where 5 s — tan - (7r/9 4) is the phase shift associ- 
ated with the potential 4 s -, vp the Fermi velocity, and 
P° = 2-kvf ^ c con duction electron density of states at 
the Fermi level. In Eq. (Q), we have also introduced the 
charge and spin boson fields, 5> c (:e) = A^($t(x) + $i(x)) 

and $ s (x) = ^($ T (a;)-$ i (x)). 

To construct soluble limits, we apply the following 
canonical transformation to the Hamiltonian, 



U 



(5) 



Using 



d$ q 



U -JT U = -JT - a2-KS(x)(2_^ aX a 



dx dx 

U + X u U = e* a *'V>X n (6) 
the transformed Hamiltonian, H e ff = U + HU, is 

H eff = H - 
J± 



Hj + AH 



H, 



Ana 



[X n FjF^e-^~ 2a ^° +H.c] 



AH= A TT - X U — 1 ) x=0 — 

np Q ax zir 



(7) 



where 5 s = ^/2S S — na. We choose an a such that what 
remains as the conduction electron part in Hj either has 
the form of a canonical fermion, or disappears entirely. 
The Toulouse points correspond to the bare values of 
the interactions such that 8 s = 0, so that the residual 
interaction AH vanishes. 
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B. Toulouse Point I of the Kondo Problem 



we can rewrite the effective Hamiltonian as follows, 



The effective Hamiltonian assumes a simple form when 
6 s is close to 6f = (3 - ^)tt, i.e., J z /4 = J}JA = 
¥ i-tan" 1 [(i - ^)tt]. Choosing a = (-y/2- l)/2 in Eq. 
(H) leads to the following transformed Hamiltnonian, 



H eff 



He 



AH 



(8) 



where 5 s = \/2(5 s — Sf) measures the deviation of the 
interaction J z from the chosen value J\. Here r\ = 
1 e -»$a (0) f j g a can onical spin- less fermion field. 

V27ra ' 

£ = X n F|F T i^, and its adjoint, d = i^F^X^, sat- 
isfy the commutation relations appropriate for a fermion, 
{d,tf} = X n +X u = 1, {d,d} = {d f ,dt} = 0. The ef- 
fective Hamiltonian, therefore, describes a resonant-level 
model pi] ]. The spin-spin correlation functions can be 
calculated using the transformed longitudinal and trans- 
verse spin operators, 



Sd 



S+ f = U + X U U 



X 



11 



(9) 



and noting that X-^ has the same dimension as 
z^F^F^". The results are 



1 



< T t S+(t)S-(0) > =< 
1 

27TO 

< T t S z (t)S z (0) > =< {d)d 

1 



1, 



:)(r) 



(dU~-)(o) > Hc 



( Po )2 
T 



V T 



(10) 



This is the usual Toulouse point pqp|,p0| , which describes 
the strong coupling, Fermi liquid fixed point. 



C. Toulouse Point II of the Kondo Problem 

We now turn to the case of S s close to J| = ^tt. 
This corresponds to an infinitely strong antiferromag- 
netic interaction, J 2 = +00. Choosing a = \/2/2 leads 
to the following terms for the transformed Hamiltonian, 
H ± = J^(X tl FlF^ + He), and AH = ^(X TT - 

*U)(^)*=05F> where = ^( SS - *!)■ Introducing a 
pseudo- fermion operator, /J, and a pseudo-boson opera- 
tor, 6^, through 



X„ 



' fafcr' 

6 f = pF^- 



(11) 



J± 



Aff-^-^K*^ (12) 



In terms of the constraint given in Eq. (J2J) can be 
rewritten as X^^o-^ = 1- At 6 s = <5|, #e// describes 
a conduction electron bath decoupled from a spin de- 
gree of freedom on which a magnetic field of strength 
h = acts along the x direction in spin space. In 

general, H^ff describes the spin-boson Hamiltonian with 

an Ohmic bath introduced in Ref . , with — being pro- 
portional to the square root of the dissipation parameter 
a defined in the macroscopic quantum coherence (MQC) 
context H|. 

Using the transformed spin operators, 



a eff 



-e^FjF^b, 



(13) 



the transverse spin-spin correlation function is given by 



<T T S+(r)S-(o) >h~(-) 2 

T 



(14) 



The calculation of the longitudinal spin-spin correlation 
function is somewhat more subtle. At 6 s = 0, there is 
only an oscillatory piece (in real time), with period 2h. 
Expanding around this point, the leading non-vanishing 
non-oscillatory term has the following long time behavior 

< T t S z (t)S z (o) > h ~ (^-) 2 (^) 2 (15) 

TTp h T 

In addition, the oscillatory piece damps out |23|] , beyond 
a time scale of ~ (jl) 2 ~\k~\- Therefore, the asymptotic 
long-time behavior of the longitudinal spin-spin correla- 
tion function has the Fermi liquid, (^r) 2 , form. This is 
the same form as that of the transverse spin-spin corre- 
lation function. Such a long-time behavior of the cor- 
relation functions is consistent with the physical picture 
that, at long times the impurity spin degrees of freedom 
is "merged" with those of the conduction electron bath. 
The (-) 2 long-time behavior for the longitudinal spin- 
spin correlation function in the spin-boson problem is 
already known in the literature |31 ,^,^3| . 

Given that this last Toulouse point occurs at an in- 
finitely strong antiferromagnetic interaction, one might 
worry about the validity of the bosonization approach. 
In Appendix B, we carry out an atomic-expansion analy- 
sis for the Kondo Hamiltonian in the limit J z >> J±, W . 
The procedure is to first diagonalize the J z coupling, tak- 
ing J± = W = 0. The lowest energy atomic configura- 
tions are 



4 



|i >= I T>d 1 4>o 
|2 | |> d | T> 



(16) 



where the subscripts d and label the impurity and the 
Wannier orbital for the conduction electrons at the ori- 
gin, respectively. J± and W couple these low energy 
configurations with other higher energy ones. Integrat- 
ing out all the high energy configurations via a canonical 
transformation leads to the following effective Hamilto- 



H eff =H' a + J^(X 12 +X 21 ) 

+J' z (Xu-X 12 )Y,<rclc lc 



(17) 



where c\ a creats a Wannier orbital of the conduction elec- 
trons at the site nearest to the origin, and J' z 



The 



same canonical transformation also leads to the following 
effective spin operators, 



1 



( S+ )eff ~ 4 T C U Xl2 



(18) 



These results from the atomic expansion, Eqs. ( |l7|) and 
(|l8|), are the direct analogs of the bosonization results, 
Eqs. @ and ©. 

To summarize, we have emphasized two aspects asso- 
ciated with the Toulouse points in the Kondo problem. 
First, we have explicitly retained the Klein factors in the 
bosonization representation of the fermion operators to 
keep track of the anticommutation relations satisfied by 
fcrmions with different spins. For the Kondo problem 
per se, the single fermion operators do not come into the 
Hamiltonian directly; only the fermion bilinear operators 
do. Therefore, identical results could have been derived 
without retaining these Klein factors. This is not true 
for the mixed valence problem, where the single fermion 
operators do appear in the hybridization term. And, as 
will be shown in detail in the following sections, it turns 
out to be essential to keep the Klein factors explicitly 
when applying the bosonization method to the mixed 
valence problem. Second, we discussed in detail a sec- 
ond Toulouse point that occurs at the infinite value for 
the longitudinal component of the Kondo exchange in- 
teraction. Through a comparison with the results of an 
expansion of the atomic limit of the Kondo Hamiltonian, 
we established that, even in this limit, bosonization can 
be applied. We also demonstrated that the correct be- 
haviors of the longitudinal spin-spin correlation functions 
at this second Toulouse point can be derived by expand- 
ing around the Toulouse point. These insights turn out 
to be quite useful to properly understand the Toulouse 
points of the mixed valence problem. 



III. THE GENERALIZED ANDERSON MODEL 

In this section we define the generalized Anderson 
model, summarize our earlier scaling results, and set up 
the bosonization formalism for the Toulouse-point analy- 
sis appropriate for the generalized Anderson model. The 
Hamiltonian of the generalized Anderson model is 

+ Y t(dU* + h.c.) + d\d a c\,c al + V p ^ c\c a 

(19) 



<7\ ,<72 ,(73 ,a"4 



Here r label the Pauli matrices, and cr =f, J.- To study 
the mixed valence regime, we focus on U — oo. The 
double occupancy configuration of the impurity is then 
excluded. The three remaining configurations, |0 > and 
|cr >= dl\0 >, have energies E Q — and E a = E%, 
respectively. The hybridization t, the density-density in- 
teraction V, and the explicit spin exchange interaction J 
describe the couplings between the impurity and the elec- 
tron bath. Anticipating the intrinsic particle-hole asym- 
metry in the mixed valence regime, we have also allowed 
for a potential scattering, V p . This Hamiltonian is gen- 
eral enough for the purpose of studying the interplay be- 
tween the local spin and charge fluctuations in the mixed 
valence problem. 

Following the procedure outlined in the previous sec- 
tion and given in more detail in Appendix A, we can 
reduce the problem to that of an impurity coupled to 
a one-dimensional non-interacting conduction electron 
bath, with one chiral component for each spin species, 
ij} a {x). For the purpose of the Coulomb gas as well as the 
Toulouse-point analyses, we also allow the longitudinal 
and transverse components of the exchange interaction, 
J z and Jj_, to take different values. The Hamiltonian can 
be rewritten as 

H = H Q + E° d X aa + H± t + Hxj + Hv 

a 

J±, v ,t, 



H 



(20) 



Here, the impurity configuration \a > runs over |0 > 
and |<7 >. This requirement amounts to the following 
constraint, 



X 



IT 



X 



ii 



X„ n — 1 



(21) 



V£ denotes the local potential that the conduction elec- 
tron of spin a experiences when the impurity is frozen 
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at the configuration \a >. In terms of the parameters of 
Eq. ©, 





V4 


J, 

4 




c= 


V- 


J, 
4 




c= 









(22) 



A. The Coulomb Gas Representation and the 
Scaling Results 

The physics of this model can be understood by focus- 
ing on the impurity degrees of freedom and tracing out 
the conduction electrons. The partition function is then 
a sum over all the histories of the impurity, each history 
being characterized by a sequence of transitions of the 
impurity configurations among |0 > and \a >. This has 
been extensively discussed in our earlier papers [jl|,|l0|,|llj . 
For completeness, we summarize the main results in this 
subsection. 

The partition function has the following Coulomb gas 
form, 



Z 
where 



E E 



exp(-S[ai,...,a n ;Ti,...T n ]) (23) 



n=0 



ai ,a n ;ti ....T n 



S[a\, ...,a n ; t\, ...t„] = y][K(cxj, a 3 ) + K(a i+ i,a J+ i) 



-K(a i} a J+1 ) - K(a i+ i,a>j)}h 



^M^ +1 ) + V^ +1 (Tt+1 , - ( 24 ) 



Here, [a%, a„; n, ...t„] labels a history, with the impu- 
rity hopping quantum mechanically from one local state, 
\ai >, to another, \ai + i >, at the (imaginary) time Tj. 
The logarithmic interaction between the hopping events 
reflects the reaction of the electron bath towards the 
changes of the impurity configurations. Here, £ is the 
ultraviolet inverse energy cutoff, and the strength of the 
logarithmic interaction is characterized by the stiffness 
constants, e t = — K(0,a) and ej — —K(a,a' ^ a) which 
in turn are determined by the bare interaction strength 
of the original Hamiltonian. Specifically, 



1 b a — r\ a r\ a — r\ a 

A 7T 7T 



(1 



where 



5% = tan- 1 (^ lC) 



(25) 



(26) 



is the phase shift that the conduction electron bath of 
spin a experiences when the impurity configuration is 
| a >. The fugacities y a ^ describe the amplitudes for a 
quantum hopping between the configurations \a > and 
|/3 >. More specifically, the charge fugacity corresponds 
the hopping amplitude between two local states with dif- 
ferent charge quantum numbers and is proportional to 
the hybridization amplitude, yo,<r = Vt — t£- Similarly, 
the spin fugacity describes the hopping amplitude be- 
tween two local states with different spin quantum num- 
bers and is given by the transverse component of the 
exchange coupling, y a ,a> = Uj = f° r o ^ a 1 . Fi- 
nally, the fields h a describe the energy splittings among 
the local configurations. In the absence of an external 
magnetic field, ho = —^E^S, and h a — \E^- 

The physical content of this Coulomb gas representa- 
tion is as follows, yj and tj are the dimensionless quanti- 
ties associated with the transverse and longitudinal cou- 
plings of the usual spin Kondo problem. y t and et, on 
the other hand, can be thought of as the analogous di- 
mensionless quantities associated with the transverse and 
longitudinal couplings in a charge Kondo problem. And 
the difference between ha- and h a , or rather the impurity 
level, E%, controls how "soft" the local charge fluctua- 
tions are; when E° d is tuned through the conduction elec- 
tron Fermi level from far below to far above, the system 
evolves from a local moment regime, through a mixed va- 
lence regime, to the empty orbital regime. In the mixed 
valence regime, Kondo-like processes in the spin channel 
and charge channel are coupled together. The partition 
function given in Eq. ( p4[ ) is the proper generalization of 
Haldane's Coulomb gas representation of the mixed va- 
lence problem such that the interplay between the spin 
and charge channels are incorporated systematically. 

Our rcnormalization group analysis of this partition 
function establishes the existence of three, and only 
three, mixed valence phases: 1) the usual strong coupling 
phase. Both the spin and charge Kondo problems are in 
the strong coupling regime; rapid fluctuations between all 
three local configurations take place and the conduction 
electrons quench both the charge and spin degrees of free- 
dom of the impurity; 2) a weak coupling phase where nei- 
ther the local charge nor the local spin degrees of freedom 
is quenched. Both the spin and charge Kondo problems 
are in the weak coupling regime, and all three atomic con- 
figurations decouple asymptotically at low energies; 3) an 
intermediate phase where the local spin degrees of free- 
dom is quenched, but the local charge degrees of freedom 
is not. Here, the charge Kondo problem is in the weak 
coupling regime despite of the fact that the spin Kondo 
problem is in the strong coupling regime. There are two 
local configurations carrying different charges which are 
decoupled asymptotically. The phase diagram is prop- 
erly given in terms of the et — £j parameter space. The 
strong coupling phase occurs when e t < 1, as well as 
over a range of tj < 1 when e t > 1- The weak cou- 
pling phase can occur only when both e t > 1 and tj > 1 . 
From Eq. (E5I) this condition means that the effective 



G 



exchange interaction has to be ferromagnetic. The weak 
coupling phase is, therefore, very likely to be of only aca- 
demic value. The intermediate phase arises over a range 
within e t > 1 and ej < 1. For the model Hamiltonian 
( |l9| ) this means an antiferromagnetic exchange interac- 
tion, and a finite attractive density-density interaction, 
V, between the impurity and co nductio n electrons. As 
will be further discussed in Sec. VIII B , taking V as an 
effective parameter this condition can be satisfied in a 
variety of realistic models. The transition between the 
different regimes is analogous to the localization phase 
transition studied in the context of MQC |23| and more 
recently in the context of transport through constrictions 
in interacting quantum wires |13[ . 



B. Construction of the Toulouse Points of the Mixed 
Valence Problem 

We now proceed to derive the Toulouse points. 
Bosonizing the conduction electrons, and making a 
canonical transformation using 

jj = e -ia$s J2 a aX <"> e - i P' p <:(I2« x <">- x °°*> e ~ i ~e i 'c (27) 
we arrive at the effective Hamiltonian, H e tt = U + HU, 
H ef f = H + E° d J2 X <™ + H t + H 3 + AH 

a 

t 



H t = 



[X ao F a e ^ ' e V2 ' + H.c. 



H,, 



J2na 

Ar^.Ftj? e -i(V5-aa)«. +Hc] 



AH=^X aa [^( d ^)' + ^&'] (28) 
^-^ irp ax 2ir irp Q ax 2ir 



2ir irp dx 



The notations are the same as in the previous section. 
&% = 5*-7raa, S c a = 5 c + 4/3- 7 ), and 5% = <%-tt(/3+ 7 ). 
Here, 5% = ^ £ CT K and 5' a = ^ E CT where 5 * is 
defined in Eq. (Eq); more specifically, 



5 s a = a(5° a - 5D/V2 

51 = {5% + 5°)/V2 

51 = (# + 5 S ) /V2 = V2S° 



(29) 



The last equality of the last equation holds in the absence 
of an external magnetic field. The latter also ensures that 
5* = 0. 

To construct a Toulouse point, we choose the param- 
eters a and (3 such that the the conduction electron op- 
erators appearing in the transformed H t and Hj have 
simple forms. The Toulouse point then corresponds to 
choosing the bare values of the interactions such that the 
transformed AH vanishes. This means taking 5i = ana, 
5% = tt(/3 + 7), and 5 C = tt( 7 - (3). From Eq. @, these 



conditions are equivalent to requiring that the phase 
shifts 5° take the following values, 

£=-L(/3 + 7 + a)7r 

^=-L(/3 + 7 -a)7r 

«5 CT =-^( 7 -/3)7r 



(30) 



According to Eq. (pa), the Coulomb gas stiffnesses at 
the Toulouse point are 



e t = J[(a + 2/3- V2) 2 + (a-2/3) 2 ] 



2a- iy 



(31) 



Eq. (|30|), together with Eqs. (|22j) and (g6|), allow us 
to determine the parameters of the generalized Anderson 
model at the Toulouse points. In terms of the phase 
shifts, they are 

5(J Z /A) ee tan- 1 ^ tan(^) - ~ tan(^)] 

= tan _1 [— tan(-7r— — — ) tan(7r ^ %= ) 

l 2 y y/2 '2 K V2 

S(V) ee tan- 1 !^ tan(^) + ~ tan(^) - tan(^)] 

, _i r l , , P + l + a 

= tan — tan 7r = 

l 2 y y/2 ' 

1 , /3 + 7 — a. 1 . 7 — /?.. 
+ — tan(7r = ) tan (n — ==-) 



5{V) = 7T 



7-/3 
V2 



(32) 



At a Toulouse point, the correlation functions can be 
calculated in terms of the single particle, off-diagonal 
spin, diagonal spin, and density operators in the trans- 
formed basis, 

(4) e// ee U^X ao U = e ia ^e 2i ^X ao 



n u - e ~" ~° x n 



s; ff = u*x n u = -- 

S: ff ^UHWaX^)U = ^aX a 



Peff 



(33) 



We found that only three independent points exist 
where the transformed Hamiltonian is exactly soluble, 
and one more point where the effective Hamiltonian as- 
sumes a simple form but is not exactly soluble. Before 
proceeding to analyze each of these points, we end this 
section with a comment on the role of the 7 term in 
the canonical transformation and, likewise, the role of 
the potential scattering term in the generalized Ander- 
son model. In general, Eq. (BOh may require that, at 
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the Toulouse points, one or several of the phase shifts, 
<5^, be outside the range specified by the unitarity limit, 
[— f , §]■ While bound states can lead to the violation of 
the unitarity bounds of the coupling constants pfj|] , in our 
case we can take advantage of the freedom of varying the 
strength of the potential scattering term to make all the 
phase shifts in Eq. ( p0| ) to fall in the range [—§,§]• This 
will be seen to be especially convenient for the purpose of 
demonstrating the consistency of the bosonization results 
with those of the atomic expansions shown in Appendix 
B. The physics is independent of the potential scattering 
term, as can be clearly seen from the 7— independence 
of the effective Hamiltonian, Eq. (^8|), as well as the 
7— independence of the Coulomb gas stiffnesses at the 
Toulouse points, Eq. (|3l|). 



IV. STRONG COUPLING TOULOUSE POINT I 

The first Toulouse point arises with the choice of a = 
& an d = At this point, Eq. @ gives e t = 
and €j = 0. The Coulomb gas analysis implies that we 
are deep in the strong coupling regime. The requirement 
that all the phase shifts 5° are within the range [—§,§] 

leads to a unique choice for 7 = — . With this choice of 
7, 8% = tt/2, 5 s a = -tt/2, and 8° = -tt/2. Equivalently, 
8(J z /4) and 8(V) are both equal to tt/2, while 8(V P ) = 
—tt/2. This corresponds to an infinite antiferromagnetic 
exchange interaction and an infinite repulsive density- 
density interaction. 

In this case, H t = --^=J2a( X ^o F cr + H.c), and 



Hj = ^(X||F/F| + H.c). We now introduce pseud- 
ofermion creation and annihilation operators /J and f a , 
and pseudoboson creation and annihilation operators b\ 
and b a , as in Eq. (|TT|). We further introduce the pseu- 
doboson creation and annihilation operators b\ and b a 
such that X ao — f^b and X ao = 6j6 Q . The effective 
Hamiltonian becomes 



H eff = 



H + E° d J2 b% + -L= + b%) 

a v a 

2np ^ 



dx 2tt 



And the constraint Eq. ( pl| ) is rewritten as 
Y J blb a + b% = \ 



(34) 



(35) 



In deriving the effective Hamiltonian, we have chosen 
5q such that the transformed potential scattering term 



vanishes. With this choice, k c and k s can be explic- 
itly written in terms of the original phase shifts as 
k c = s/2(SZ + 8%) and n s = ^2{8° - 8% - tt). They 
vanish at the Toulouse point. 

At the Toulouse point, the Hamiltonian describes a 
three-level system decoupled from the conduction elec- 
trons. It is the three-level analog of the two level spin- 
boson Hamiltonian with zero dissipation. The three-level 
system can be diagonalized exactly, with three eigen- 
states, |s0 >= u\A > +v\o >, \t >= \B >, and 
\sl >= —v\A > +u\o >. Here, 

\A>=-^(b\+b\)\vac> 

\B >= -L(6{ - b\)\vac > 
\o>=b\\vac> (36) 



define a new basis for the atomic states. In the absence 
of hybridization, these three configurations have ener- 
gies E A = E°- ij, E B = E° d + and E = 0. 
The hybridization term mixes \A > and |o >, lead- 
ing to the bonding and anti-bonding states |s0 > and 



\sl >. 
Id" 



u and v are the coherence factors, u 2 

E^ — J^/Aira 



1 



). The energies of the 



y<(E°-J ± /4:Tra) 2 +8(t/V2^) 2 

eigenstates are E s0 = [{E% - J^/Attq) 2 + 8(i/v / 2^) 2 ]/2, 
E t = E° + J ± /4ttci, and E sl = [(E° d - J ± /4ttci) + 

\J{E° d - J^/Atto) 2 + 8(t/V2^) 2 ]/2. Given that J z is 
antiferromagnetic, Jj_ should also be antiferromagnetic 
(positive). Therefore, irrespective of E d , |s0 > is always 
the ground state. 

The ground state of the whole system is \gs >= 
U + \s0 > \FS >. The correlation functions can be cal- 
culated from the transformed single particle, spin and 
density operators 

(4W/ - e^/^e^/^ 2 ^=(X Ao + aX Bo )F$ 

S^ff = e lv/ ^* s -{X AA — X BB + X BA - X AB ) 



1 



Peff — (X AA + X BB — X ao ) 



(37) 



The single-particle and the transverse spin susceptibility 
are straightforward to calculate and have the Fermi liquid 
forms, 



< T T( f CT (T)4(0) > ~ £2 

T 

< T t S-(t)S+(0) > ~ (^) 2 

T 



(38) 



At the Toulouse point, the longitudinal spin-spin correla- 
tion function and the density-density correlation function 
again have the oscillatory behaviors in real time at the 
Toulouse point. As for the longitudinal spin correlation 
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function near the Toulouse point II of the Kondo prob- 
lem, a direct expansion in terms of the deviation from 
the Toulouse point leads to the following Fermi liquid 
behaviors for these two correlation functions, 



clear that both the exchange interaction (J±) term and 
the hybridization (t) term favor the same singlet state 
\s0 >= u\A > +v\o >, for arbitrary values of E d . No 
level crossing occurs. 



< T t S z (t)S z (0) > ~ ( 



2-Kpohs T 



<T rP (T)p(0)>~(-^-) 2 (^) 2 (39) 
2irp h c t 

where h s — J±/4ira and h c = tj\phm. This Toulouse 
point, therefore, describes a strong coupling, Fermi liquid 
state. 

Unlike for the Kondo problem, keeping track of the 
anticommutation relation between fermions of different 
spins in the boson representation plays an essential role. 
Failing to do that, the sign of Hj term would be reversed, 
leading to the atomic configuration \A > having energy 
E°, + 4^- instead of E°, — -rS and \B > having energy 

E°, — 4^- instead of E°, + -p*-. In that case, the spin- 
et 47ra a 47ra 1 " 

flip exchange interaction (Jj_) would make the config- 
uration \B > energetically more favorable than \A >, 
while the hybridization (t) term, which mixes \A > with 
|o >, would favor |s0 > instead of \t >= \B >. This 
competition between the spin-exchange and hybridiza- 
tion would then lead to a level-crossing as E° d is varied: 
the ground state changes from \t > to |s0 > as the d— level 
varies from far below to far above the Fermi level. Thc 
level-crossing found in earlier works on related problems 
^8| , ^9| results from this failing to keep track of the an- 
ticommutation relation between the bosonized fermion 
fields. 

We can gain much physical insights into our results 
by determining the quantum numbers of the respective 
impurity eigenstates. Within thc bosonization approach, 
the meaning of the atomic configurations, \A >, \B > 
and |o >, is somewhat obscure. The physical content of 
these configurations becomes transparent once we com- 
pare them with the atomic configurations that appear 
in a perturbation expansion of the original Hamiltonian 
in terms of J±/J z , J±/V, t/J z , t/V, W/J z , and W/V. 
This atomic expansion is a natural procedure given that 
Jz,V ^> J±,t,W at the Toulouse point. The details of 
this atomic expa nsio n can be found in Appendix B. From 
Eqs. ( p6|) and (p7|), \A >, \B > and \o > are identified 
with 



\ A> = ^(\ t>d I !>o -I l>d I T>o) 
\B> = ^=(\ T>d I l>o +| l>d I T>o) 
|O>H0> d |2> 



(40) 



\A > is the local singlet formed between the impurity 
spin and the conduction electron spin at the impurity 
site, \B > the S z = state of the local triplet, and |o > 
the singlet with the impurity empty of electrons and the 
local conduction electron orbital doubly occupied. It is 



V. STRONG COUPLING TOULOUSE POINT II 



We now choose a = 4p and (3 



V2-2 



The cor- 



responding Coulomb gas stiffnesses are e t = 1/2 and 
tj = 0, and the Coulomb gas analysis would again pre- 
dict that the system is in the strong coupling phase. 
As in the previous section, the requirement that all the 
phase shifts fall in the range [— f , f ] leads to a unique 

choice for 7 = | — 4 ■ This choice corresponds to 

5% = tt/2, 51 = -tt/2, and 8% = ^^tt. Equivalently, 

6(J z /4) - it/2, 5{V) = -^P±tt, and 8{V p ) - ^±tt. 
We have an infinite antiferromagnetic exchange interac- 
tion and a large, but finite, attractive density-density 
interaction. 



In this case, H t = -£= Y.AXaoF^e 



H 



3 4ira 
t = i?t 1 



H.c), and 

(X^iFTF^ + H.c). Introducing a fermion field, 
*<=, and defining /t = X a0 F a F^, we derive 



thc following effective Hamiltonian 



Hf fl = E kVtVk + E° d flf a + t[Q2 ft)n + H.c] 



(41) 



It is convenient to introduce a new basis for the atomic 
configurations, \A >= ^(/j + /])|0 >, \B >= ^(/{ - 

/|)|0 >, and |o >. Thc corresponding atomic levels are 

E a = E° d -{^ a ,E B = E° + fa, and E = 0. In this 
new basis, the constraint becomes 



Xaa + Xbb + X 



1 



(42) 



(43) 



and thc effective Hamiltonian can be rewritten as 
H eff = J2 EkvW + V2t(X Ao?7 + H.c.) 

k 

+ [E° - ^]X AA + [E° d + ^]X BB 

Among the three impurity configurations \A >, \B >, 
and |o >, > and |o > hybridize with the conduc- 
tion electrons, while \B > is decoupled. Therefore, this 
Hamiltonian can be diagonalized exactly. Thc Hilbcrt 
space factorizes into two sectors. Again, since J z is anti- 
ferromagnetic, Ji should also be taken as antiferromag- 
netic. It is clear that, the low lying excitations lie in the 
sector in which the ground state is 



Cf) = \o > (f> + \A > (j) A 



(44) 







where <fi a and cj) A are the conduction electron wave func- 
tions such that 4> is the ground state of a resonant level 
model with X Ao treated as a free fermion operator. The 
ground state energy (relative to a free Fermi sea) is 
that of the resonant-level model with a resonance width 
2itp t 2 and an effective d— level E° d — 4±-. 

Using 4 = e iV2* s <r/2 e i(V2-2)^/2 X ^ g+ = 

e iVJ*.jf Uj S* = {X AB + X BA )/2, and p = (X AA + 
Xbb — X do ), and noticing that X Ao has the dimension of 
2 «0c wc g nc ] thg^ , A \\ ^jjg correlation functions have 



the Fermi liquid form, 



< T T d a (r)4(0) > 

7 

<T t ^-(t)5 + (0)>-(^) 2 



< T t S z (t)S z (0) > ~ ( 



K s \2rPo\2 
2TT Po h s T 



< T T p(r)p(0) > ~ (^) 2 

T 



(45) 



where k s is again the deviation from the Toulouse point. 
This establishes the strong coupling, Fermi liquid nature 
of this Toulouse point. 

Except for the change of sign in J±, the effective 
Hamiltonian ( j4l| ) is identical to that of Refs. and p9| . 
The sign change for Jj_, however, plays a crucial role. It 
implies the absence of a level-crossing as a function of E3. 
Similar to what happens in the previous subsection, the 
result is the Fermi liquid behavior for the longitudinal 
spin-spin correlation function, as for the transverse spin- 
spin correlation function. No unexpected fixed points as 
discussed in Refs. and occur. The bosonization 
results are consistent with our previous renormalization 
group results. 



A TOULOUSE POINT FOR THE 
INTERMEDIATE PHASE 



VI. 



In Refs. |9|,|10|,y_l[ we found a new mixed valence phase, 
which we called the intermediate phase. It is a non-Fermi 
liquid with quasiparticlc-like spin excitations and inco- 
herent charge excitations. In this section, we present a 
Toulouse point which clearly exhibits the physics of the 
intermediate phase. It occurs with the choice of a = & 

and P = — The corresponding Coulomb gas stiff- 
nesses are tt = 1 and e 7 = 0. While it is not possible to 
determine the precise boundary between the intermediate 
phase and the strong coupling phase, it is not inconsis- 
tent with the Coulomb gas results that these values of 
the Coulomb gas stiffnesses lie close to such a boundary. 
Taking 7 = specifies the phase shifts at the Toulouse 
point, 8% = § , 8% = -§ and 8° 



8{J Z ) = f , 8{V) = -f , and 8{V P ) 
ters correspond to an infinite antiferromagnetic exchange 



■| . Equivalently, 
^. These parame- 



interaction and an infinite attractive density-density in- 
teraction. 

In the canonically-transformed Hamiltonian, H t — 



/27ra *-i<tL~""- " - ' j J -*3 47raV^ 1 U- t l 

H.c). The effective Hamiltonian can be written as 
H eff =J2 EkC ka c ka + V2tV2^[X Ao c^c L + H.c] 



ka 



J±_ 

4ira J 



)X 



AA 



J± 
Ana 



BB 



+ -^—{X AA + Xbb - X 00 )(cjc T + cjej.) 

^Po 

+ ^(X AB + X BA ){c\ C] - c\ Cl ) (46) 

Po 



Here, 



A >= jsZA-cX'oFt) and \B 



>= 



^Efft - ^"^!)' Again, with the requirement on 8q 
such that the transformed potential scattering term van- 
ishes, k c and k s can be written explicitly as k c = 
v2(<$£ + (5^)and n s = \2{8° — 8° — n), and are non-zero 
only away from the Toulouse point. 

In this effective Hamiltonian, the charge sector is de- 
scribed by a genuine charge Kondo model. \A > and 
\o > play the role of | f> and | |> of the spin Kondo 
problem and should be thought of as objects carrying 
charge 2 and 0, respectively. The transformed hybridiza- 
tion term is the direct analog of the spin-flip term in the 
spin Kondo problem, with a change of the charge quan- 
tum number by two replacing a change of the spin quan- 
tum number by 1 in the latter. The residual interaction 



in the charge sector, 



2irp 



is the analog of the longitu- 



dinal exchange interaction in the spin Kondo problem, 
with the density playing the role of the spin in the lat- 
ter. In the conventional notation, ~ j^ har s c = 2iy / 7ra and 
1 jchargc _ Ke__ rj^ esS ential difference between the 

4 % 27Tp a 

charge Kondo problem in this mixed valence context and 
the spin Kondo problem lies in the symmetry-breaking 
field. In the latter, the spin symmetry guarantees that 
no explicit magnetic field term will be generated in the 
absence of an applied field. In our charge Kondo prob- 
lem, the particle-hole symmetry is explicitly broken, and 
the symmetry-breaking field h chavgc = E° d — is in gen- 
eral non-zero. For the impurity problem, the condition 
that the renormalized /i char s° vanishes can be achieved 
only through fine-tuning the bare d— level E° d to a criti- 
cal value E c d . 

When h cha - lgc = is enforced, a zero temperature 
quantum phase transition takes place as k c is increased 
through zero. The transition is characterized by a 
Kosterlitz-Thouless transition in the charge sector; the 
spin sector is not critical. The phenomenology of the 
intermediate phase is recovered on the negative k c side, 
to which the remaining of this section is devoted. Here, 
the charge sector is described by the weak coupling fixed 
point of the charge-Kondo problem, while the spin ex- 
citations by the strong coupling, Fermi liquid-like fixed 
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point of the Kondo problem. A spin-charge separation 
takes place. 

Within the charge sector, the impurity configuration 
in the ground state is entirely \o > for h chars ° < 0, and 
\A > for h chm ' gc > 0. This is the result of infinite charge 
susceptibility in the corresponding ferromagnetic charge 
Kondo problem. Exactly at h cha - Tgc — 0, namely, when 
E° is tuned to the critical value E c d = j^, the impurity 
degrees of freedom in the ground state involve an equal, 
incoherent, mixture of |o > and \A >. Schematically, the 
ground state wavefunction can be written as 

<j> = \A > <t> A + \o >(j> (47) 

where <Pa and <f> are the wave functions of the conduction 
electrons such that <\> solve a ferromagnetic Kondo model 
with zero magnetic field. With /j char g° = 0, the interme- 
diate mixed valence dynamics applies at all temperatures. 
When E° d is moved away from the critical value, a finite 
cross-over temperature T co ~ \E d — E d \ emerges. The in- 
termediate mixed valence dynamics continue to apply at 
T > T co . At low temperatures (T < T co ), however, the 
charge fluctuations become gapped out. Such a cross- 
over can already be inferred from the renormalization 
group trajectories in our previous work p0| . 

The correlation functions that describe the critical 
dynamics associated with the intermediate mixed va- 
lence phase can be calculated explicitly. Consider first 
the single particle Green function. Using (dj^e// = 
gis./^e-i*./^^)^ + X Bo )F«, we find that 

< T T d CT (r)4(0) > ~ (^)[*+5(i-^) 2 ] (48) 

T 

The fact that the exponent is interaction (/t c )-dependent 
has already been anticipated by the Coulomb gas analysis 
|p|, ^0|JTl1l . What is new here is the explicit demonstration 
that the exponent of the single particle Green's function 
in the intermediate phase is the sum of two terms, one 
from spin excitations and the other from charge excita- 
tions. The spin contribution is |, independent of inter- 
actions. The charge contribution, on the other hand, is 
|(1 — \/2^f) 2 and varies as the interaction strength is 
changed. The exponent is always larger than one, sig- 
naling a quasiparticle residue vanishing in a power-law 
fashion in the intermediate phase - a non-Fermi liquid 
indeed. 

For the spin-spin correlation functions, using = 
F\F ie ^^{X AA - X BB + X AB - Xba), and S z eff = 
^(Xab+Xba), and following a procedure similar to that 
used in calculating the spin-spin correlation functions in 
the second Toulouse point of the Kondo problem, we de- 
rive the following, 

<T t S-(t)S+(o) >~ (^) 2 

T 

< T t S*(t)S z (o) > ~ (^^) 2 (^) 2 (49) 
2irp h s t 



where h s = -p^. This establishes the Fermi- liquid be- 
havior of the spin-spin correlation functions. 

We now turn to the density-density correlation func- 
tion, using p e ff = Xaa + Xbb — X 00 . Due to the 
crgodicity breaking, there is a disconnected contribu- 
tion to the density-density correlation function that is 
t— independent. The connected piece is algebraic, which 
we found to be, 

< T tP (t) P (o) > conncctod ~ ^) (^)(-4kc) (50) 

This algebraic piece, once again, has an interaction- 
dependent exponent. We note that, the fact that the 
exponent vanishes with k c has already been noted before 
in related problems |34|,|35| . 

Other correlation functions in the charge sector also 
have an algebraic behavior with interaction-dependent 
exponents. For instance, the excitonic correlation func- 
tion, using (4c CT ) e// = e- i ^"(^)(X Ao + aX Bo )F^F l , 
has the following form, 

<T r (^ C td CT )(r)(^4c CT )(0)>^(^) 2(1 -* )2 (51) 

a a 

Finally, we consider the pairing susceptibility. Follow- 
ing a similar procedure, we find that 

< T T (^c ff 4)(r)£44)(0) >~ (^) ( ^ )2 (52) 

a a 

which is enhanced compared to the Fermi-liquid case. 
This makes it plausible that the ground state in the corre- 
sponding lattice model is superconducting. In that case, 
the intermediate mixed valence dynamics would describe 
the physics in the normal state, i.e., at temperatures be- 
tween the transition temperature and some upper cutoff 
energy scale. 

To summarize, the intermediate phase has spin-charge 
separation, a quasiparticle residue vanishing in a power- 
law fashion, self-similar local correlation functions with 
interaction-dependent exponents. These characteristics 
bear strong similarity to those of the Luttinger liquid in 
one dimensional interaction fermion systems p5| , |26[ |. 

VII. A POINT WITH SIMPLE, BUT NOT 
EXACTLY SOLUBLE, EFFECTIVE 
HAMILTONIAN 

The three Toulouse points that we have discussed 
for the mixed valence problem is exactly soluble. In 
this section, we discuss one more point in the inter- 
action parameter space which is described by a sim- 
ple effective Hamiltonian. It arises from the choice of 
a = — (1 — and j3 — i^j^- The Coulomb gas stiff- 
nesses are £t = \ and Cj = 2. In this case, there is a 



11 



range of 7, 



x/2 
4 



< 7 < 



iV2 



1. that satisfies the re- 



quirement that all the phase shifts 8^ fall in the range 

and (Jq = (— I + ^)tt- This range corresponds to a large 
ferromagnetic J z and a large repulsive V (and attractive 



V^). For instance, choosing 7 = y„ 



= — ^ corresponds 



to *(J,/4) = 
chossing 7 = 7 

5(F) = f , and 5(V P ) = - 

Independent 
the choice of 7, Ht — - 



^±tt, S(V) = f , and 5(V p ) = - f , while 



1 leads to S(J z /4) = -§, 



V2-1, 



- »<i> 



irrelevant, we can take the Jj_ = 



As if. 



of 

- ff-c), 
is strongly 



limit and consider 
Hj as a perturbation later. Introducing a fermion field 



»7 = ^: 



'2iva 



and 



>>= Icr > F„F%, the effective 



Hamiltonian has the following form, 

H eff = E + Wo + X °-^ + H - C l - E d X °° ( 53 ) 



This Hamiltonian has the same form as that investi- 
gated analytically and numerically in our earlier studies 
of a spinless two band model in large dimensions f{6|| . 
It has also been found pTj] to arise within a particu- 
lar rotation scheme in the bosonized form of the Ander- 
son impurity model with additional screening channels. 
What distinguishes this effective Hamiltonian from those 
of the previous subsections is the fact that, the operator 
(-Xfo + X [) is now a rank-2 matrix. In fact, H^tt is 
a rank-2 generalization of the Emery-Kivelson resonant- 
level model for the Toulouse point of the two-channel 
spin— i Kondo problem ^0|. The non-abelian nature of 
the phase space makes the problem non-quadratic and 
not exactly soluble. In the following, we establish the 
nature of the solutions in the two limits, —E° d >> \t\, W, 
aadE%»\t\,W. 

Consider first the limit of —E% >> \t\,W. Here, the 
configurations | t> and l> lie at low energies, while |o > 
should be treated as a high energy configuration. The hy- 
bridization t term mixes the low energy and high energy 
configurations, and can be eliminated with a canonical 
transformation. This results in an effective Hamiltonian, 
Ho + (-T&r)(*tT - XlOtfv - 1/2). The sign of the 

1 d ' 

effective interaction makes the spin-flipping term even 
less relevant. The result is a ground state with a double 
degeneracy. 

For the opposite limit, E° d >> \t\, W, \o > alone lies at 
low energies, while | f> and J.> are high energy configu- 
rations that can be eliminated. The resulting low energy 
behavior is a potential scattering problem. The ground 
state is a singlet, and the low lying excitations have the 
Fermi liquid form. Therefore as the impurity level moves 
from far below to far above the Fermi level, the system 
evolves from a non-Fermi liquid with a doublet ground 
state to an empty-orbital Fermi liquid. 



We can in fact understand the doublet character of 
the ground state in the limit of — E° d » \t\, W already 
from the renormalization group equations for the origi- 
nal Hamiltonian, Eq. (|l9|). The relevant renormalization 
group equations are given in the Appendix B of Ref. ( 
|10|). For a large negative E d , the contribution to the 
scaling of the spin-flip Ji term from the second order 
in hybridization t term is small. The scaling of Ji is 
then entirely determined by J 2 , as in the Kondo prob- 
lem. Given that J z is ferromagnetic, the system remains 
a spin doublet at the fixed point. As E° d moves closer to 
the Fermi level, the contribution of the second order in t 
term to the scaling of J± becomes more important, and 
eventually dominates the J z contribution. These conclu- 
sions are consistent with the numerical calculations we 
have performed on the Hamiltonian (|53|). 

Yet one more line of reasoning leads us to the same 
physical picture of the phases of the Hamiltonian (p3). 
We have shown earlier |36| that the Hamiltonian ( |53| ) 
maps onto the spinless resonant level model with particle- 
hole symmetry and in the limit of vanishing W/V r \ and 
W/t T \ (where W, V r \ and t t \ are, respectively, the conduc- 
tion electron bandwidth, the interaction strength, and 
the hybridization of the resonant level model). In this 
mapping, the parameter t of Eq. ( |55| ) is of the order of 
W, while E d = V r i/2 — t r i. For a fixed t r i, increasing E d 
in the Hamiltonian (|5^) is equivalent to increasing V r i 
in the resonant level model with particle-hole symmetry. 
When t r i/W is small, the resonant level model under- 
goes a phase transition from a doublet to a singlet as V r i 
increases from a strongly attractive value to a large repul- 
sive one. The mixed valence regime of the Hamiltonian 
( |53| ) is the strong coupling version of the critical regime 
associated with this phase transition of the resonant level 
model. Whether the physics of the resonant level model 
at large values of t r i /W and V r i /W is smoothly connected 
to the physics of the resonant level model at small t r i/W 
and Vri/W, or a phase transition separates these two 
regimes, remain unclear and is a problem deserving fur- 
ther investigation. 



VIII. CONCLUDING REMARKS 
A. Related Works 

The mixed valence problem is closely related to prob- 
lems in the macroscopic quantum tunneling (MQT) and 
macroscopic quantum coherence (MQC) where a few lo- 
cal degrees of freedom are coupled to a bath of low energy 
excitations p^ ]. In the MQT and MQC literature, typ- 
ically only two-level systems are considered. The mixed 
valence problem amounts to a three-level generalization. 

The mixed valence problem has historically been stud- 
ied in the context of impurities embedded in a metal. It 
was stressed early on that, in order to satisfy the Friedel 
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sum rule, screening interactions between the impurity 
and additional channels of conduction electrons have to 
be introduced p3|. Similar screening interactions have 
also been studied in the contexts of two-level systems 
and macroscopic quantum tunneling [39|j30| ] . 

In the context of high temperature superconductivity, 
it has been suggested [Q that the mixed valence regime 
of an impurity model, with a repulsive impurity density- 
conduction electron density interaction and a large num- 
ber of screening conduction electron orbitals, would be a 
local model exhibiting the novel non-Fermi liquid char- 
acteristics of the marginal Fermi liquid, i. e. logarithmi- 
cally divergent charge and spin susceptibilities. The nu- 
merical results of Rcf . |l7j indicate divergent charge and 
spin susceptibilities near the mixed valence point. How- 
ever, the numerical results are not conclusive because the 
susceptibility enhancement may be just due to a crossover 
to a local moment regime as the impurity level is varied. 
The divergence is especially difficult to see in the spin 
susceptibility, as the susceptibility continues to increase 
when E° d is decreased through the transition regime [jl7| . 

Within the renormalization group analysis, the effect 
of screening channels is to modify the initial conditions 
of the renormalization group flow [pj. The additional 
screening channels are passive observers which slow down 
the response of the impurity degrees of freedom (hence 
increase the orthogonality) but do not participate in the 
formation of fixed points other than those within the 
renormalization group classification. This is consistent 
with the considerations of Giamarchi et. al. (4^] and 
Guinea et. al. pC( |. From this perspective, the Anderson 
model with additional screening channels has the same 
low energy behavior as the generalized Anderson model 
considered here as well as in our earlier works p, ^0| , pT| . 



Recently, however, two groups [|18 19 have sought to 
study an exactly soluble point of the Anderson model 
with additional screening channels. These works, using 
the bosonization method, have reached conclusions that 
would signal a novel fixed point unexpected from our 
previous renormalization group results. One of the ma- 
jor conclusions in this work - the arguments for which 
being detailed in Sees. [j"v| and [v| - is that the fixed 
points discussed in Refs. llsUTij are the result of using 
incorrect bosonization expressions for the fermion oper- 
ators. What was missing are the Klein factors that keep 
track of the anticommutation relations of the fermions 
of different spin species. When the correct bosoniza- 
tion expressions are used, the results become compatible 
with our previous renormalization group results. We are 
therefore forced to the conclusion that, the only known 
generic solutions to the spin— i mixed valence problem 
in the generalized Anderson impurity models are the 
three phases identified within our previous renormaliza- 
tion group scheme. 



B. The Intermediate Phase in Realistic Models 

The intermediate mixed valence state that we have 
identified represents a new state of strongly correlated 
electron systems. It will therefore be interesting to ex- 
plore the phenomenological consequences of this non- 
Fermi liquid state, and to search for such a state in other 
strongly correlated electron models. In this connection, 
two important issues need to be addressed. 

First, within the generalized Anderson model that we 
have studied, the interaction parameter regime for this 
intermediate phase corresponds to a range of antifer- 
romagnetic exchange interaction, and finite attractive 
density-density interaction, between the impurity and the 
conduction electrons [in addition to the large repulsive 
on-site Hubbard interaction]. It is important to address 
how this attractive density-density interaction between 
the impurity and conduction electrons can be generated 
from purely repulsive interactions in more realistic im- 
purity and lattice models. Several directions have been 
explored along this direction. Additional screening chan- 
nels enhance the orthogonality effects (io) and effectively 
play the role of attractive density-density interactions. 
Alternatively, dynamical screening effects that arise from 
integrating out high energy configurations in strongly 
correlated electron systems can generate effective attrac- 
tive density-density interactions From the perspec- 
tive of a local approach such as in the limit of infinite 
dimensions, this effect is especially important since the 
effective bandwidth of the self-consistent conduction elec- 
tron bath is usually quite small |ll|,E^]. The generation 
of attractive density-density interaction from purely re- 
pulsive interactions have also been discussed in a class of 
multi-band models Q . 

Second, in general non-Fermi liquid states in impurity 
models can be realized only through fine-tuning parame- 
ters. For the intermediate mixed valence phase, there is 
one parameter that needs to be fine-tuned, namely the 
impurity level E° d . This places constraints on the realiza- 
tion of the intermediate state in real materials displaying 
impurity physics. However, the situation becomes much 
better in real materials that are described by lattice mod- 
els, as alluded to at the beginning of this paper. In lattice 
models, the mixed valence condition can be satisfied over 
a range of densities H,[HJ . 



C. Summary 

In this paper, we have studied several exactly solu- 
ble points in the mixed valence regime of the general- 
ized Anderson model. We found three such points and 
clarified the physics of each. Two of these points have 
the strong coupling, Fermi liquid behavior. The third 
Toulouse point describes the novel intermediate mixed 
valence phase. Our explicit results about the correlation 
functions at this last Toulouse point clarified the exci- 
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tation spectra in this new phase of strongly correlated 
electron systems. Finally, we established that, once the 
anticommutation relations between the bosonized forms 
of the fermion fields are taken care of, the Toulouse-point 
results derived within the bosonization formalism become 
consistent with our previous renormalization group re- 
sults. 

From a more general perspective, our results clearly il- 
lustrate the importance of the competition between local 
charge and spin fluctuations in strongly correlated elec- 
tron systems. In the intermediate phase, the competing 
low energy local charge and spin fluctuations lead to a 
spin-charge separation, Fermi-liquid-like spin susceptibil- 
ities, non- Fermi- liquid charge susceptibilities, non- Fermi 
liquid single particle spectral functions, and an enhanced 
pairing susceptibility. This represents a new route to- 
wards spin-charge separation alternative to what leads 
to the Luttinger liquid in the interacting fermion models 
in one dimension. In our case, the underlying physics is 
local in nature, while in one dimension, it is dominated 
by long wavelength fluctuations. Nonetheless, the char- 
acteristics of the correlation functions in the intermediate 
phase have strong similarities to those of the Luttinger 
liquid. 
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APPENDIX A: BOSONIZATION FOR IMPURITY 
PROBLEMS 

In this appendix, we summarize the bosonization pro- 
cedure relevant to our discussion. For the purpose 
of studying an impurity problem with contact interac- 
tion at r = only, the S— wave component of the con- 
duction electrons alone needs to be kept. This S— wave 
component can in turn be written as a superposition of 
an incoming component and an outgoing one defined on 
the half-axis r G [0, +oo), with the boundary condition 
that these two components be equal at the origin pi) . 
This boundary condition allows us to keep only one, say 
the incoming, component while simultaneously extend- 
ing the problem to the full axis, x G (-co, +oo). Denot- 
ing this component as tp a [x), discarding the uncoupled 
higher-than-S 1 — wave components of the conduction elec- 
trons, and linearizing the energy dispersion about the 



Fermi level, we can rewrite the Hamiltonian for the non- 
interacting conduction electron bath, H , as follows, 



Vp 



dxilpa — — 

ax 



(Al) 



where vf is the Fermi velocity. The linearized conduction 
electron dispersion assumes the form Ek = vp(k — kp). 
The density of states is p(e) — J2k &{Eh — e) 



27TVf 



We can now introduce a boson representation of the 
■0-field. To keep track of the anticommutation relation 
between the fermion fields, we write the boson represen- 
tation of the fermion fields as follows, 



i 



(A2) 



Here a is a cutoff scale which can be taken as a lat- 
tice spacing. $(x) is defined in terms of the Tomonaga 
bosons, 



. . 2irx , T 



E 

q>0 



2n ( -hi 
+ib„„t 



iqx—qa/2 



-iqx—qa/2\ 



(A3) 



In Eq. (A3), b qa and b qa are the Tomonaga bosons, 




k+q a^ko 



k 



(A4) 



Notice that § a (0) involves only the q ^ components of 
the Tomonaga bosons. L is the length of the dimension: 
x G [-L/2, L/2]. N a and F a represent the zero modes 
of the boson fields. N a is the deviation of the conduc- 
tion electron occupation number from the ground state 
value. It represents the q = counterpart of the finite q 
Tomonaga boson occupation number n qa = b qa .b qa . The 
boson Hilbert space is spanned by \N a ,n qa >. Within 
this Hilbert space, the operator FJ raises N a by one, 



while its adjoint, F a , lowers N a by one ]2^ 
are traditionally called the "Klein factors" 
the following relations, 

F*F„ = F a Fl = 1 
F*F 9 = -F s Fl 
F F- — — F-F 



T 



These 
ley satisfy 



(A5) 



They commute with b q<J and b qtT , for q ^ 0. 

In this boson rep resentation, the non-interacting 
Hamiltonian Eq. ( |Al|) has the following form, 

f+OO 



and the density operator 

p a (x) = iplip a = 



dx 



1 rf$«r 

27T dx 



(A6) 



(A7) 
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APPENDIX B: ATOMIC EXPANSION 

Consider first the To ulous e point II of the Kondo prob- 
lem discussed in Sec. II C . Given that the vicinity of 
this Toulouse point corresponds to J z » J±, W, we can 
carry out an expansion in terms of 4^ and W-. For this 
purpose, we rewrite the original Kondo Hamiltonian as 



N 



H = wY] (cU +lCT + H.c.) + -±(S + c\ C] + sA 



(Bl) 



where cL creates a Wannier orbital at site i in the Wilson 



basis [|2|. To the leading order, we need simply to diag- 
onalize an atomic problem defined in the Hilbert space 
of the impurity spin doublet and cJ CT . The lowest energy 
atomic states are given 



|i>=IT>dU>o 

|2>=||> d |T>o 



(B2) 



For finite J± and W, these low energy atomic states 
become coupled with the high energy atomic states, 
\Ao >= \a > d |2 > , \Bo >= \a > d |0 > , and 
\Ca >— \a >d \<J >o- Integrating out these high energy 
states via a canonical transformation 



W 

X 



+X Al . 2 c n + X Bl . 2 cL - H.c.) 



leads to an effective Hamiltonian -He// — e s He . 
the second order in W/ J z , 

H e 'ff = H'o + J±-{Xi2 + ^21) 



To 



where 



+j' z (Xn - X22) } J <rc\ a c 1(T 



N 

2 — 1, IT 



(B3) 



(B4) 



ii'- 

J, 



and J' z 

The effective spin operators are , S^j — and 
are given by 



S. 



eff 
tit 



- ~ s S z e- s 



1 



(Xn ^ X22) 



St,= e s S + e 



S = 2(^)V 1T c u X 12 

'J Z 



(B5) 



Eqs. ( |F33j - |B5D are the direct analogs of Eqs. ( |12[ ) and 



We now turn to the strong coupling Toulouse point I 
of the mixed valence problem discussed in Sec. IV. The 
low energy atomic states are, 



|i >= I T> d 1 4>o 

|2>=||> d |T>o 
\0>= |0> d |2> 



(B6) 



The mixed valence condition amounts to the requirement 
on E d such that the energy difference between |1 >, |2 > 
and \0 >, 2E' d = E° d + V + \V P \ - i, is small. The 
canonically transformed Hamiltonian is 



H. 



Ai 

eff 



H'o + E d (Xn + X22 — Xoo) + Ja_{X\2 + X21) 
+t{X 10 -X20 + H.C.) 

+J' z (Xn-X22)J2' JC l^ 



+ V^X n +X22)^2 C l^lc 



(B7) 



where V' z 
are 



The canonically-transformed operators 



+ WW* 

(d\)eff ~ (— + y)c[,X 10 

+ w w + 



V 



a eff 



7 eff 



1 



(Xn — X22) 



Peff ~ (-^11 + -^22 — Xqo) 



(B8) 



These expressions are the direct analogs of the expres- 
sions derived using the bosonization representation, Eqs. 
©and©. 

Finally, we consider the intermediate Toulouse point 
of Sec. VI. The low energy atomic states are, 

\A >= u-^= °W >d \o >o +v|0 > d |2 > 
\B >= 7/| X)k >d |S >o 

\0 >= |0 > d |0 > (B9) 

The canonically transformed Hamiltonian is 

a 

eff 



H°L =H' a + J ± [-X AA + X BB ] + t'{X A0 c l]Cll + H.c) 



+V / (X AA -X o)CE c Lcic 



(BIO) 



where t' 



(l + V 2 )W 2 rrn, • 11 

^-^ — '- — 1 he canonically 



transformed S* f f , Sfyy, and p e // are similar to those 
given in Eq. (B8), while 
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1 w w 

(4) e// -=( — + —) Cls (aX AO + X BO ) (Bll) 
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